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Abstract 

We study the heat kernel of the sub-Laplacian L on the CR sphere S 2n+1 . An 
explicit and geometrically meaningful formula for the heat kernel is obtained. As 
a by-product we recover in a simple way the Green function of the conformal sub- 
Laplacian — L + n 2 that was obtained by Geller [12) . and also get an explicit formula 
for the sub-Riemannian distance. The key point is to work in a set of coordinates that 
reflects the symmetries coming from the fibration S 2n+1 — > CP™. 
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1 Introduction 



The purpose of this work is to study the heat kernel of the sub-Laplacian of the stan- 
dard CR structure on § 2n+1 . More precisely, we will be interested in explicit analytic 
representations and small time asymptotics of the kernel. 

A key point in our study is to take advantage of the radial symmetries of the fibration 
g2n+i _^ Qpn £ introduce coordinates that are adapted to the geometry of the problem. 
In particular, it is shown that the kernel has a cylindric invariance property and actually 
only depends on two variables r, 9. The variable 9 is the local fiber coordinate of the 
fibration and r is a radial coordinate on CP™ . In these coordinates the cylindric part of 
the sub-Laplacian L is the operator 

d 2 d d 2 

L = —r + ((2n - 1) cot r - tanr)— + tan 2 r— -. 
or z or o9 z 

Eigenvalues and eigenvectors for L are computed and as a consequence of the general 
Minakshisundaram-Pleijel expansion theorem we deduce that the subelliptic heat kernel 
(issued from the north pole) can be written 

k=— oo m=0 

where Pm is a Jacobi polynomial and X m) k = ^m{m + \k\ + n) + 2|A;|n. This formula 
is very useful to study the long-time behavior of the heat kernel but seems difficult to 
use in the study of small-time asymptotics or for the purpose of proving upper and lower 
bounds. In order to derive small-time asymptotics of the kernel, we give another analytic 
expression for pt(r,0) which is much more geometrically meaningful. This formula is 
obtained thanks to the observation that the Reeb vector field T of the CR structure of 
g2n+i commu t es with the sub-Laplacian. This commutation implies that the subelliptic 
heat semigroup e tL can be written e~ tT2 e t ^, where A is the Laplace-Beltrami operator of 
the standard Riemannian structure on § 2?l+1 . The formula we obtain is explicit enough 
to recover Geller's formula (see [12]) for the fundamental solution of the conformal sub- 
Laplacian — L + n 2 . By using the steepest descent method, it also allows to derive the 
small-time asymptotics of the heat kernel. A by-product of this small-time asymptotics is 
a previously unknown explicit formula for the sub-Riemannian distance. 

To put things in perspective let us observe that the study of explicit expressions for 
subelliptic heat kernels has generated a great amount of work (see [l], [2] [4], [3], [5] 
[7], [11] and the references therein). The motivations for finding explicit formulas are 
numerous, among them we can cite: sharp constant in functional inequalities (see [6], 
|14j). computation of the sub-Riemannian metric (see [3]), sharp upper and lower bounds 
for the heat kernel (see [5], [TO]), and semigroup sub-commutations (see [13]). However, 
despite such numerous works, very few explicit and tractable formulas are actually known 
and most of them are restricted to a Lie group framework. The present work gives explicit 
and tractable expressions that hold in a natural sub-Riemannian model. 
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2 The sub-Laplacian on S 2n+1 

2.1 Geometry of the standard CR sphere 

We consider the odd dimensional sphere 

§2n+l ={z= ( Zu ... f Zn+l) £ C n+ \ \\Z\\ = 1}. 

It is a strictly pseudo convex CR manifold (see [9]) whose geometry can be described as 
follows. There is a natural group action of S 1 on § 2n+1 which is defined by 

Ol,--- ,z n ) (e l6 zi,--- ,e ie z n ). 

The generator of this action shall be denoted by T throughout the paper. We have for all 
/ e C°°(S 2n+1 ) 

d 

~d~9' 

so that 



Zjf dzj Z ^ dzj 



This action induces a fibration (circle bundle) from S 2n+1 to the projective complex space 
CP n . The vector field T is the Reeb vector field (characteristic direction) of the pseudo- 
Hermitian contact form 



. n+l 

r] = --^{zjdzj - zjdzj). 

3=1 



2.2 The sub-Laplacian 

For j = 1, • • • , n + 1, let us denote 

where S = Y^k=i Zfe af^> an< ^ define the second order differential operator L on C°°(S 2n+1 ) 
as follows: 

n+l 

i=i 

It can be checked that L is the CR sub-Laplacian of the previously described structure 
(see for instance [6], [8]). It is essentially self-adjoint on C co (§ 2n+1 ) with respect to the 
uniform measure of S 2n+1 and related to the Laplace-Beltrami operator A of the standard 
Riemannian structure on §> 2n+1 by the formula: 

L = A — T 2 . 
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We can observe, a fact which will be important for us, that L and T commute, that is, on 
smooth functions TL = LT. 

To study L we now introduce a set of coordinates that takes into account the sym- 
metries of the fibration S 2n+1 — > CF n . Let (wi,--- ,w n ,6) be local coordinates for 
S 2n+1 , where (w\,--- ,w n ) are the local inhomogeneous coordinates for CP n given by 
Wj = Zj/z n+ i, and 6 is the local fiber coordinate, i.e., (wi,-- - ,w n ) parametrizes the 
complex lines passing through the north poleS, while 6 determines a point on the line that 
is of unit distance from the north pole. More explicitly, these coordinates are given by 

(„,..., Wn , „) I - 7 _, . . . , - 7 _, I , (2.2) 



where p = y^j=i l^il 2 ; ^ ^ M/27rZ, and w £ CP n . In these coordinates, it is clear that 
T = -j^. Our goal is now to compute the sub-Laplacian L. In the sequel we denote 

n r-. 

Proposition 2.1 In t/ie coordinates we have 

t=*a+ p 2 ) E ^ + 4(1 + p2)R ^ + " 2 5F - 2 * 2 + " *>£■ 

fc=i 



Proof. From (|2.ip . we know that 



By using now the diffeomorphism 

(>!,••• ,W n ,9,K) - 



nO if) i 



and then restrict to the sphere on which we have: 



9 

K=l,— =0, 
OK 



we compute that on § 2n+1 , for 1 < k < 



■n 



9 ^-5 c -id 9 



9zfc ' dw k 

9 = y/TT?* 9 -^- 



1 We call north pole the point with complex coordinates Zi = 0, • • • , z n +i = 1, it is therefore the point 
with real coordinates (0, • • • ,0, 1, 0). 
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and 



Tedious but straightforward computations lead then to 

k,j=l 3 

which implies that 

n+l 



^ \ ZjZk dz.m + ZjZk &z-dz k ) - 2 89 2 ■ 

;,7=1 J J 



Moreover, it is not hard to compute 

n+1 / d 2 d 2 \ ( d 2 — 1 d 2 1 — 9 



fc=i 



where 7£ = Ylk=i w k~ef^- Finally, it is easy to see that 5 + 5 = 0. Hence we have the 
conclusion. □ 

Remark 2.2 Notice that T = j^, and thus the Laplace- Beltrami operator is given by: 

A r 9 2 

W n 

= ^ + p 2 )iZo^d=+ ^ + + a + p 2 ) w 2 - 2 ^ 2 + w - 

k=l 

Due to the symmetries of the fibration S 2ra+1 — > CP™, in the study of the heat kernel, 
it will be enough to compute the radial part of L with respect to the cylindrical variables 
(P,0). 

Let us consider the following second order differential operator 

~ , 2X 2 d 2 f(2n- l)(l + p 2 ) , 2x \ d 2 d 2 



which is defined on the space T> of smooth functions / : R>o x R/27rZ — > R that satisfies 

dp 



— — if p = 0. It is seen that L is essentially self-adjoint on V with respect to the 



„2n-l 



measure — e T dpd6. 

( l +p2) 2n+l 
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Proposition 2.3 Let us denote by ip the map from § 2n+1 to M> x M/2vrZ such that 
For every f £ T>, we have 

L(fo^) = (Lf)o^. 
Proof. Notice that by symmetries of the fibration, we have 

(sr 92 ) ff o vo = ( (-— + 2n ~ 1 9 

diVkdvoj: ) \\Adp 2 4p dp 



\k=l 

and 

n(f o v) = n(f ^) = (Q^) /) ° ^ 

Together with Proposition 12.11 we have the conclusion. □ 

Finally, instead of p, it will be expedient to introduce the variable r which is defined 
by p = tanr. It is then easy to see that we can write L as 

d 2 d d 2 

L= — J + ((2n-l)cotr-tanr)— + tan 2 r— -. (2.3) 
or 1 or oO z 

This is the expression of L which shall be the most convenient for us and that is going to 
be used throughout the paper. 

We can observe that L is symmetric with respect to the measure 

27T n 

du r = (sin r) 2 " - " 1 cos rdrdO. 

r(n) 

The normalization is chosen in such a way that 

2vr n+1 



^ [ 2 dp r =p(§ 2n+1 ) 
J -it JO 



nJo v ' T(n + 1) 

Remark 2.4 In the case of § 3 (n = 1), which is isomorphic to the Lie group SU(2), we 
obtain 

r m 2,2 d 2 ({1 + P 2 ?\ d 2 8 2 
L = (l + p 2 ) 2 ^+[- — hr+P 



dp 2 \ p J dp d9 2 

d 2 d 2 d 2 

— T + 2cot2r— + tan r— ? 
dr z or ot) z 



This coincides with the result in J3J/. 
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Remark 2.5 By Remark \2.2\ we see that the radial part of the Laplace- Beltrami operator 
in cylindrical coordinates is 

d 2 did 2 
A r = — + ((2ra - l)cotr - tanr)— + 



dr 2 dr cos 2 r d9 2 

On the other hand, since in the coordinates \2.2\) . we have z n +\ = cosre ie , it is clear that 
the Riemannian distance form the north pole 5, satisfies 

cos 5 = cos r cos 9. 

An easy calculation shows that by making the change of variable cos 5 = cos r cos 9, the 
operator A r acts on functions depending only on 5 as 

° 2 „ ? d 
W + 2nC ° U d5- 

This expression is known to indeed be the expression of the radial part of A in spherical 
coordinates. 



3 The subelliptic heat kernel on § 2n+1 
3.1 Spectral decomposition of the heat kernel 

From the expression of L above, it is not hard to see that the kernel of Pt = e tL issued 
from the north pole only comes from the radial part L and depends on (r, 6). We denote 
it by p t (r,9). 

Proposition 3.1 Fort > 0, r £ [0, 9 E [— 7r,7r], the subelliptic kernel has the following 
spectral decomposition: 

k=— oom=0 

where X m ,k = ^m(m + \ k\ + n) + 2|/c|n and 



pn-l,|fc|f \ _ V ^' (f-i _ T \n-l+m(-t , \\k\+m\ 

m K) ~2™ m \{\-xy-\\ + x)\K\dx™ KK } 
is a Jacobi polynomial. 

Proof. The idea is to expand pt(r, 9) as a Fourier series in 9. Let 

+oo 



Pt(r,9)= ^ e <fce fc (t,r), 



fc=— oo 
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be this Fourier expansion. Since pt satisfies -Jr = Lp t , we have 

d(f> k d 2 cf) k d(f> k 2 2 

— — = - n + ( zre — 1 ) cot r — tan r ) — k tan rcp^.. 

By writing (p k (t,r) in the form 

fc (t,r)=e- 2 "l fe l t (cosr)l fc lr 7fc (t,cos2r), 

we get 



where 



= 4* fc fe), (3.4) 

tf fe = (1 - x 2 ) + [(|fc| + 1 - n) - (|fc| + 1 + n)x]|-. 

In fact (|3.4p is well-known as the Jacobi differential equation, and the eigenvectors are 
given by 

(_-i\m jm 

P"- l >\ k \(x) = — ^—((1 - x) n - 1+m (l + xp +m ), 

which satisfies that 

^ k (PZ- m )(x) = -m(m + n+ ^P^^ix). 
Therefore, we obtain the spectral decomposition 

+oo +oo 

p t (r,9) = ^a mifc e-[ 4m ( m+ l fc l +n ) +2 l fc l n ]*e ifce (cosr)l /£ lp™- 1 'l fc l(cos2r), 

k=— oo m=0 

where the constants a mtk 's are to be determined by the initial condition at time t = 0. 

To compute them, we use the fact that (P^T m (x)(l + x)l fc l/ 2 ) m > is an orthogonal 
basis of L 2 ([-l, !],(!- x) n ~ l dx), i.e., 



/l 
P^- 1 'l fc l(x)P"" 1 ' |fc| (x)(l-x)"- 1 (l+^) |fc| ^ 



2 n+\k\ T(m + n)T(m+ \k\ + 1) 

"ml- 



2m +\k\+n T(m + l)T(m + n+\k\) 
For a smooth function f(r,8), we can write 

+oo +oo 

/M)= E E 6 ^ fcejP ™" 1,|fc| ( cos2r )-( 1 + cos2r ) |fc|/2 

fc=— oo m=0 

where the {6fc )m }'s are constants, and thus 

+oo +oo 



/(o,o)= y: e 6 ^" 1 ' 1 ^ 1 ) 2 '" 172 - 



k=— oo m=0 
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Now, since 

K 

Pt (r,9)J^e)d f , r 

Pt(r, 0)f(r, 0)(sinr) 2n_1 cos rdrdO 

-r^J TOO „IL 

EE/ / a m ,^ mifc 2l fc l/ 2 e- A -^|P^- 1 'l fe l| 2 (cosr) 2 l fc l+ 1 (sinr) 2 "- 1 dr^ 

._ n J -IT JO 



■K JO 

2ir n r7T 



T(n) 



7T Jo 



- „ +00 +OO nnr „ — 

T(n) 

k=—oo m=0 ' 

9 „n +°° +°° / r- 



E ^aM.tWe^^-""!'!/ 2 - 1 y I y 2 n+ l fc l +1 (co S r) 2 l fe l +1 |i^- 1 'l fc l| 2 ( S inr) 2n - 1 dr J d0 



T(n) 

fc=— oom=0 



E E a m , fc W^ W 2- ri - |fc|/2 - 1 (2 7 r)||Pr 1 ' |fc| H 2 



T(n) 

k=— oom=0 

where A m) fc = 4m(m + |fc| + n) + 2|&|n, we obtain that 

lim r [ 2 Pt fdfi r = f (0,0) 

t ^ J-TT JO 

as soon as a m , k = ^(2m + | k\ + n) . □ 

The spectral decomposition of the heat kernel is explicit and useful but is not really 
geometrically meaningful. We shall now study another representation of the kernel which 
is more geometrically meaningful and which will turn out to be much more convenient 
when dealing with the small-time asymptotics problem. The key idea is to observe that 
since A and ^ commutes, by Remark 12.21 we formally have 

e tL = e- 4 ^e* A . (3.5) 

This gives a way to express the sub-Riemannian heat kernel in terms of the Riemannian 
one. Let us recall that the Riemannian heat kernel writes 

V( \ + °° 

% (cos<5) = E(™ + ^- m{m+2n)t Ocos5), (3-6) 

m=0 

where, as above, 5 is the Riemannian distance from the north pole and 

C»M = (~ir r(m + 2n)r(n+l/2) 1 d™ 

m{ ' 2™ T{2n)T{m + l)r(n + m + l/2)(l-x 2 ) n - 1 / 2 dx m{ ' 

is a Gegenbauer polynomial. Another expression of qt(cos5) which is useful for the com- 
putation of small-time asymptotics is 
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1 (J — 2fc7T)^ 

where V{t, 5) = ^= JJ fcgZ e « . 

Using the commutation (|3.5p and the formula cos (5 = cos r cos we then infer the 
following proposition. 

Proposition 3.2 For t > 0, r £ [0,vr/2), 9 € [-n,ir], 

1 /'+ 00 (y+ifl) 2 
Pt(r, 9) = e a qJcosr coshyjdy. (3.8) 

V4vrt J-oo 

Proof. Let 



e 44 ft (cos r cosh y)dy, 



oo 
d 



and Lq = -T^z + ((2n — 1) cot r — tanr) J-p, then we have 



L = Ln + tan r 



90 2 ' 



Using the fact that 



d_ 
at 



V 

e 4t 



'Ant 



dy 2 



v 

e « 



'Ant 



and 



we get 



0_ 

Of 



Lh t (r, 9) = (L + tan 2 r-^ ) h t (r, 9) 



(ft (cos r cos 9)) = A(ft(cos r cos 9)), 



2 



r+oo 






/ (y+iS) 2 \ 
e n | 




L 


i< 


1 —oo 




i VAwi J 



+oo 



+oo 



A 



e n 



i d 2 

cos 2 r a# 2 
dq t 



d 2 

ft(cosrcoshy) | + tan 2 r-^2 



e « 



'Ant 



qt (cos r cosh y) 



_(y±wy 
e n 



d 2 

ft(cosrcoshy) ) + tan 2 r-^j 



_(y±wy 
e it 



dy 



qt (cos r cosh y) 



dy 



+ 



1 die n 



(y+jey 



I Ant j dt cos 2 r eft 



2 d I e it 
q t - tan r— 
at 



4vrt 



ft 



On the other hand, it suffices to check the initial condition for functions of the form 
f(r, 9) = e tX9 g(r) where A£l and g is smooth. We observe that 

h t (r,9)f(r,9)d^ r = e tx \e t ^g)(0). 



□ 



o J o 



Thus ht(r,9) is the desired subelliptic heat kernel. 
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Proposition 3.3 For A G C, ReA >0, r£ [0,tt/2), (9 G [-ir,n], 

+ ~ „<„ = r , — -im ^ (3 . 9) 

J -co 2"+2 7r n+i C osh s/y 2 + 4A - cos r cos(6» + iy) 



Proof. Since 

f + OO /"+OC 



\/47T 7- 

We want to compute 

n + OO 2 



p t (r,e)e- n " t -7di= / e-T Pt (r,6)e- n2t dt 
o Vo 

/ / e « q t [cosr cos{& + iy))—=dy 

J-oo Jo yt 



r°° -nH-£+^ t r n • w dt 

/ e « q t (cosr cos{u + iy))—=.. 

Jo yt 



Notice that 



Vi Jo y/i 

c - > A 2 +4Ay / ^A 

Vn 2 - A 

However, by the result in Taylor [[IS], pp. 95], we have that 

A~ x e~ tA = ^7r-(" +1 )r(n)(2cosht - 2 cos ,5)-™, 
where A = \Jn 2 — A and cos (5 = cos r cos(# + iy). Plug in f = y/y 2 + 4A, we obtain 

1 -V^+4AV^A rfn) 



V" 2 - A 2 n + 1 7r™+ 1 (cosh yV + 4A - cos r cos(6> + iy) 



hence complete the proof. □ 
We can deduce the Green function of — L + n 2 immediately from the above proposition. 

Proposition 3.4 The Green function of the conformal sub-Laplacian —L + n 2 is given by 
G{r,9) - [ ' 2 ' 



8vr n+1 (1 - 2 cos r cos 9 + cos 2 r) n l 2 
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Proof. Let us assume r ^ 0, 9 7^ 0, and let A — ^ in (j3.9|) . we have 



G(r,0) 



r(n) /-+ 00 ciy 



2n+2 7r n+i (cosh y — cos r cos(# + iy)) T 

T(n) r+°° dy 



2«+2 7r n+i y ^ (cosh y(l — cos r cos 6) — i cos r sin 9 sinh y) T 
T(n) 1 /'+ 00 dy 



2 n + 2 7r n+1 (1-2 cos r cos 6* + cos 2 r)™/ 2 7_ 00 (cosh y) 



Notice that 



+oo 



dy vr(n-2)!! 



, (cosh y) n (n — 1)!! 

2(n-2)!! 
(n- 1)!! 

where n!! denotes the double factorial such that 



when re > is odd 
when n > is even 



re!! = re • (re — 2) • • • 5 • 3 • 1 when re > is odd 
= re • (re — 2) • • • 6 • 4 • 2 when re > is even 



Moreover, since 



( n — 2)n 2 n_1 /re\ 2 
rfnH = r - whenre>0isodd 



(re-1)!! vr V2 



re\ 2 



_ 2 n 2 r ^— j when re > is even 
we obtain 

/•+°° // ? / /r?\2 

This implies our conclusion. □ 
Remark 3.5 This agrees with the result by Geller in llStf . 

3.2 Asymptotics of the subelliptic heat kernel in small times 

First, we study the asymptotics of the subelliptic heat kernel when t — > on the cut-locus 
of 0. From (|3.7p . we already know the following asymptotics of the heat kernel 

1 / d \ _sl f ( 2 n ( n ~ l)(sin<5 — (5 cos 5) \ 2\ 



q t (cos5) = e -4t 1+ re 2 -^ >\ . '- )t + 0(t 

(4nt) n+ 2 \sm6J \ \ 5 2 sm5 J j 

(3.10) 

where 5 S [0, n). Here 5 is the Riemannian distance. Together with (|3.8p . we first deduce 
the following small-time-asymptotics of the subelliptic heat kernel on the diagonal. 
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Proposition 3.6 When t — > 0, 



1 



/ " ,,UJ:i "= (^TT^« + jB " t + ^)' 



where A n = /_« -4—,/,, and S n = (l + (n 2 + ^^Sg^ )) dy. 

Proof. We know that 

1 f°° _y^_ 

p t (0,0) = —== e ^ qt (cosh y)dy. 

Plug in (|3.10p . we have the desired small time asymptotics. □ 
Proposition 3.7 For 9 £ (0,vr) ; t ->■ 0, 



on— 1 



^ ( °^ ) = 23^(n-l) ! e "^ (1 + 0(t)) 



Proof. Let 6* G (0, 7r), we have 



1 

Pt{0,6) = / e 4t % (cosh(y - z0))dy 



Mvrt 



By Cauchy integral theorem, this is the same as integrating along the horizontal line in 
the complex plane by shifting up i6 from the real axis, i.e., 

1 f°° (y+i6) 2 

p t {Q,6) = —= I e a q t (coshy)dy 

V47Tt J-oo 

Moreover, by (|3.10p . we know that 



This gives 



g t (coshy) ~ t ^ tt • , e 4t 

(4TTt) n+ 2 Vsmhy/ 



ell f°° y n _vye 



By the residue theorem, we get 

(iyS \ 
e -2ty n 
- — , — kni 
{smhy) n 



I (sinhy) n 
1 d n ~ l 

-27TZ > -, — — T 

^ (n-l)\dy n - 1 



e -^2 n y n (y + kTTiy 
(ev - e~y) n 



y=—k7vi 
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Write W(y) 



(y+kni) 



(ev-e-v) 



fyz, W{y) is analytic around —km, and satisfies 



W{-km) 



^ -^kn<2n ' 



Hence the residue is 



1 & 



n-l r 



(n- 1)! dy n ~ l L 



e-^2 n y n W(y) 



y=—kni 



iyB 



This is a product of e at and a polynomial of degree n — 1 in 1/t. We are only interested in 
the leading term which plays the dominant role when t — > 0. Thus we have the equivalence 



-2vri d 



n—l 



(n - 1)! dy^ 1 
At the end, we conclude 

Pt(o,e) 



e-^2 n y n W{y) 



kn+n n+ 1 nn—1 



y=—kni 



^-5-0 



{-l) kn+n n 



[n 



1 )l 2 n-2 t n-l 



kir8 

e 2t 



e « 



(_ 1 Nfcn+n 7r n+l 6 ,n-l 

> - - e 2t 



that is 



gn— 1 



Pt(0,f) 



2 3n t 2n( n _ 



e — (1 + 0(t)) 



□ 



Come to the points that do not lie on the cut-locus, i.e., r / 0. First we deduce the 
case for (r, 0) . 

Proposition 3.8 For r G (0, we have 

Pt(r, 0) 



e 4t 



(4vrt) r 



n +5 V sin TV V 1 — rcotr 



1 



os t — y 0. 

Proof. By proposition 



Pt(r,0) 



together with (|3.1U|) . it gives that 

Pt(r,0) ~t^o 



(4vrt) 



e 4t q t (cosr coshy)dy, 



(Ji(t) + J 2 (t)) 



n+l 
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where 

/ ■/) , ( 

'coshj/<-J— \ y/l — COS 2 V COsh 2 y 



f ;< 2 + (arccos(cosrcoBha)) 2 / aiCCOsfcOS T COSh if) \ 

J H*) = / . e 4t I n I ^ 

J CO; 



and 

. , f V 2 -(cosh" 1 (cos r cosh y)) 2 ( COSh" 1 (COS T COSh V ) \ 

Mt) = / e « I dy 

Jcoshj/>^^ y y cos r cosh y—lj 

We can analyze Ji(£) and J2W by Laplace method. First, notice that in [— cosh - ^^—), cosh _1 (^^:)] , 

f(y) = U 2 + (arccos(cosr coshy)) 2 
has a unique minimum at y = 0, where 

/"(0) = 2(1 - rcotr). 
Hence by Laplace method, we can easily obtain that 



Vsmr/ VI 



Ant 



r cot r 

On the other hand, on (—00, — cosh _1 (— -!— )) U (cosh -1 (— |— ), +00) , the function 

V * COS 1 ' / \ * COS / / 

g{y) = y 2 ~ (cosh -1 (cos r cosh y)) 2 

has no minimum, which implies that J2(t) is negligible with respect to Ji(t) in small t. 
Hence the conclusion. □ 

We can now extend the result to the 9^0 case by applying the steepest descent 
method. 

Lemma 3.9 For r G (0, §), 9 £ [-vr, vr], 

f(y) = (y + *#) 2 + (arccos(cos r cosh y)) 2 

defined on the strip |Re(y)| < cosh -1 (^7) feos a critical point at iip(r,9), where <p(r,9) 
is the unique solution in [—it, it] to the equation 

. „. . m arccos(cos(/9(r, 6>)cosr) 

<p( r ,0) + = cosrsiny(r,g) \ , i - 

V 1 — cos^ r cos z </?(r, 0) 

Proof. Let u = cos r cos 99, 

9 / . arccos(cos 09 cos r) \ sin 2 r / u(r, 9) arccos u(r, 9) \ 

— 09 - cos r sin u> , = = ; — —7 1 1 — = — - 

v / l-cos 2 rcos 2 a9y l-u(r,0) 2 y ^l-^M) / 

is positive, thus 9 = <p — cos r sin 99 arccos ( cos v? c °s>-) j s a bijection from [— vr,7r] onto itself, 



-cos 2 r cos 2 <p 



hence the uniqueness. □ 
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Moreover, observe that at (p(r,9), 



.it,. , ... 2sin 2 r / u(r, 9) arccos u(r, 
f W(r,e)) = - j—^z 1 



l-u(r,9) 2 y ^l-n 2 (r,^) 

is positive, where u(r,9) = cos r cos ip(r, 0). 

By using the steepest descent method we can deduce 

Proposition 3.10 Let r 6 (0, §), 6 G [— 7T, tt]. Then when t ->■ 0, 

1 (arccos u(r, #)) n e 4tsin2 M r . e )) 



Pt(r,0) 



J t->o 



(4-7ri) n+ 2 sinr A u(r,e) arccos u(r,fl) (i _ u ( r . g) 2 ) U 2 



where u(r,9) = cos r cos <^(r, #). 

Remark 3.11 By symmetry, the sub-Riemannian distance from the north pole to any 
point on § 2n+1 only depends on r and 9. If we denote it by d(r,8), then from the previous 
propositions, 

(1) For 8 € [-vr,vr], 

d 2 (0,9) = 2n\9\ - 9 2 

( 2 ) For 9^ [-7T, it], r £ (0,f), 

d 2^ e) _ ^{r,9) + 9)HB^r 
sin 2 (ip(r, 9)) 

In particular, the sub-Riemannian diameter of § 2n+1 is it. For a study of the sub- 
Riemannian geodesies on §> 2n+1 ; we refer to fTt 
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